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Nambu-Jona-Lasinio (NJL) model with Kobayashi-Maskawa-'t Hooft (KMT) term is one of low 
energy effective theory of QCD which includes the Ua(1) anomaly. We investigate nonet meson 
properties in this model with three flavors of quarks. We employ two type of regularizations the 
dimensional and sharp cutoff ones. The model parameters are fixed phenomenologically for each 
regularization. Evaluating the kaon decay constant, the r\ meson mass and the topological suscepti- 
bility, we show the regularization dependence of the results and discuss the applicability of the NJL 
model. 
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I. INTRODUCTION 

The unique QCD Lagrangian is found by imposing the 
local SU C (3) gauge symmetry, Lorentz invariance, local- 
ity and renormalizability in four space-time dimensions. 
It is believed to be a fundamental theory of strong inter- 
action between quarks and gluons. One of the features 
of QCD dynamics is "asymptotic freedom" which legiti- 
mates perturbation theory at short distances. Unfortu- 
nately non-perturbative effects cannot be avoided in the 
confinement phase which takes place at low energy scale, 
i.e. where the QCD coupling is not small. 

Nambu-Jona-Lasinio (NJL) model [l| is a well-known 
and often used low energy effective theory of QCD [2|- 
4]. Nambu and Jona-Lasinio have introduced a four- 
fermion interaction to describe the attractive force be- 
tween fermions. In this model, chiral symmetry is spon- 
taneously broken by non-vanishing expectation value for 
a composite operator constructed by the fermion and 
anti-fermion fields and the fermion mass is dynamically 
generated. The NJL model and its generalizations are ex- 
tremely useful in the study of the light meson properties 
at low energy. 

The four-fermion interaction is a dimension six oper- 
ator in four space-time dimensions therefore the model 
is non-renormalizable and depends on the regularization 
procedure. In order to regularize fermion loop integrals 
one usually introduces a momentum scale A to cutoff 
integration momenta higher than A. Another regulariza- 
tion, the dimensional one, is an analytic regularization; 



one calculates loop integrals as analytic functions of the 
space-time dimensions and uses them for the value of the 
space-time dimensions less than four JEHU [l2[ • 

In the present paper we study NJL model in the dimen- 
sional regularization and compare the results with ones 
obtained in the sharp cutoff regularization. NJL model in 
both regularizations describes well the dynamical break- 
ing of chiral symmetry and the tt and a meson properties 
in vacuum [8[. Yet at finite density (especially when the 
cutoff scale is close to the Fermi momentum) important 
contributions to the quark loop integrals are dropped in 
the cutoff regularization scheme @, while dimensional 
regularization leads to results consistent with QCD even 
in the region of asymptotic freedom, therefore we find 
a strong regularization dependence in the extended NJL 
model with an attractive force in the color anti-triplet 
channel for a large chemical potential [l2j]. Such a high 
density state may be realized in astrophysical objects. It 
is expected that the model can be tested by observing 
the structure of dense stars. 

Non-negligible regularization dependence is observed 
even when considering the system at T = fi = 0. The 
cutoff scale is usually taken to be lower than rf meson 
mass (ttw — 958MeV) to reproduce light meson proper- 
ties 0, 0, EH • Then the m,,/ in the low energy effective 
theory is not well-defined quantity. A is phenomenologi- 
cally fixed and it is considered as a scale above which the 
effective model may lose its validity. Furthermore, the 
cutoff regularization may break some symmetry of the 
Lagrangian. The cutoff regularization may cause some 
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unexpected effects in the system where the strange quark 
and E/a(1) anomaly play an important role. Therefore we 
launch a plan to study the three flavor system with Ua(1) 
anomaly by using the dimensional regularization. 

In the present paper the extended NJL model including 
the Kobayashi-Maskawa-'t Hooft (KMT) term (HCDl is 
regarded as a low energy effective theory of QCD with 
Ua(1) anomaly. We consider three-flavor light quarks, 
u, d and s, and investigate the nonet meson proper- 
ties. In Sec. II. we introduce the model Lagrangian and 
briefly review regularization schemes. In Sec. Ill we cal- 
culate the meson masses, decay constants and topologi- 
cal susceptibility in the leading order of 1/N C expansion. 
The results depend on the regularization parameters. In 
Sec. IV we evaluate the kaon decay constant, the 77 meson 
mass, rrirj, and the topological susceptibility ^asa func- 
tion of the space-time dimension D associated with loop 
integrals. In Sec. V the dependence of physical quanti- 
ties on the current up-quark mass, m u , is studied. After 
phenomenologically fixing the model parameters, we dis- 
cuss the two regularizations dependence of the results. 
This paper shows that the model Lagrangian is not sat- 
isfactory in either of regularizations. Some concluding 
remarks are given in Sec. VI. 



II. NJL MODEL WITH U A (l) ANOMALY 

The NJL model is one of the simplest models to de- 
scribe the dynamical symmetry breaking. It is often 
used to study the light meson properties at low en- 
ergy scale. Since the ?7a(1) anomaly induces the mass 
difference between r\ and rj mesons, the NJL model 
should be extended to include the contribution from the 
C/a(1) anomaly in evaluating the nonet meson system 
constructed out of the three flavor light quarks. 



A. Model set up 

Kobayashi and Maskawa have introduced an interac- 
tion written in the determinant of the composite operator 
constructed by quark and anti-quark fields in the nonet 
meson system [l5j. In QCD #-term can be transformed 
to the determinant term by ?7a(1) transformation. Thus 
we can include the contribution of the J7a(1) anomaly 
through the determinant term in the NJL model. In the 
present paper we consider the four- and six-fermion in- 
teraction invariant under SUl(3) <8> SUr(3) global flavor 
symmetry and start with the Lagrangian, 



where 



C 



NJL 



Qj + £4 + £-6 



(1) 



^2 H^aQj 



a=0 

£ 6 = -K [det&(l - 75)% + h.c. ] , 



^2 H^Kqj 



(2) 
(3) 



the subscripts are the flavor indices, to expresses the 
current quark mass matrix, A a are the Gell-Mann matri- 
ces in the flavor space, G and K represent the effective 
coupling constants for four- and six-fermion interaction, 
respectively. The determinant in Cq concerns the matrix 
elements labeled by the flavor indices, ij. The order of 
the coupling constants are supposed to be GN C ~ O(l), 
KN% ~ O(l) therefore we work in the framework of the 
1 /N c expansion. In this paper we do not care about the 
flavor mixing and set the mass matrix to have a diagonal 
form, to — diag(m„, nid, m s ). The current quark mass 
explicitly breaks the global SU(3) flavor symmetry. Be- 
low we consider the SU (2) isospin symmetric case and 
take m u = for simplicity. 

The chiral condensates (uu) , (dd) and (ss) generate the 
constituent quark masses, m* , m* d and to* inside mesons. 
To evaluate the constituent quark mass we solve the gap 
equations which are derived by differentiating the ther- 
modynamic potential with respect to to*, m* d and to*. 
The gap equation is obtained in the leading order of 1 /N c 
expansion 

m* u = m* d = m u + AG{i tr S«) + 2K(z teS d )(i tr 5 s ), (4) 
m* =m s +4G(itrS s ) + 2K(itrS u )(itrS d ), (5) 

s v s,c 1 v s,c 1 K s,c 1 v y 

where the symbol tr stands for the trace in spinor and 
color indices. trS 11 represent the chiral condensates, 
— itrS" 1 = (uu) and — itr5 s = (ss), which are given by 
the trace of the quark propagator inside mesons, 



-itrS 1 = N r tr 



s 1 (p) 



d D p 
i(2ir) L 



s 1 (p), 



(6) 



1 



— to* + ie 

Here we indicate the space-time dimension for internal 
quark fields by D. Bel 
notational simplicity. 



quark fields by D. Below we omit the subscripts in tr for 



B. Regularization schemes 

The fermion loop integral in Eq.([6|) is divergent in four 
space-time dimensions to obtain a finite result we have 
to regularize it. The four- and six-fermion interaction 
are written in terms of the dimension six and nine oper- 
ators, respectively. Thus the operators in £4 and Cq are 
irrelevant in four dimensions. It means that the results 
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depend on regularization procedures. Here we use two 
different procedures; one is the three-momentum sharp 
cutoff regularization and the other is the dimensional reg- 
ularization. 

In the three-momentum sharp cutoff method, we cut 
off the space momentum component integral above the 
scale, A, 



JD 



P 



dpo 
2tt 



1 d 3 p 
(2tt) 3 



(7) 



In the dimensional regularization scheme, we regularize 
the divergent integral with the help of analytic continu- 
ation of the integral as a function of the space-time di- 
mension D to a non-integer value less than four, 



id 



P 



2(4tt) 



-D/2 



(27T) 



D 



T{D /2) 



dpp 



D-l 



(8) 



This is a kind of an analytic regularization. We regard 
the space-time dimensions D in the fermion loop integral 
as one of the parameters of the effective model of QCD. 
The dimensional regularization is applied to momentum 
integrals only for internal fermion lines. 

The common parameters of the models considered here 
are the coupling constants G and K, the current quark 
masses m u (= m^), m s . In the cutoff scheme, the cut- 
off scale A is one more parameter. On the other hand 
in the dimensional regularization, we consider the space- 
time dimension D as one of the model parameters. In 
this case we have to introduce one more parameter, the 
renormalization scale M , to obtain results with the cor- 
rect mass dimension. Thus the parameters in these two 
regularization methods are aligned as follows, 

Cutoff: G, K, m u (= m^), m s , A, 
Dimensional: G, K, m u (= rrid), m S) D, Mq. 

All the parameters should be fixed phenomenologically. 



III. MESON MASS AND DECAY CONSTANT 

In this section we shall evaluate the properties of the 
nonet meson system. Here we calculate the meson mass, 
meson decay constant and topological susceptibility in 
the two regularizations. 



where the index a denotes the isospin channel and P 
stands for the meson species, ir and K. The flavor- 
dependent effective couplings K a are defined by 

K 3 = G+ ^KitrS 3 , for tt°, (10) 

K 6 = G+^KitrS u , for K°,K°. (11) 

In the leading order of the l/N c expansion the self-energy 
for each meson, P, is given by 



n P (fc 2 )<5 a/ 



d D P 

i(2n) L 



tr[ 75 T a S> + k/2) l5 Tj i S 1 ( P - fc/2)] , 



(12) 



where the trace runs over flavor, spinor and color indices. 
The SU(3) matrices, T a , corresponding to different chan- 
nels are, T 3 = A 3 for tt°, T e = (A 6 + iX 7 )/V2 for K° and 
t\ for K°. Thus the self-energy for tt°, K° and K° read 

TL n (k 2 =ml) = 2nr(fc 2 -m 2 ), (13) 
Il K (k 2 = m 2 K ) = 2IVf{k 2 = m 2 K ), (14) 

where iTf (k 2 ) is the loop integral 

If(fc 2 ) - J-^tr[ l5 S\p + k/2) l5 S3(p-k/2)] 
1 / i trS 1 i trS j \ 



2 \ m* 



(15) 



with 



d D p 



trl 



/y(fc2) J i(2n) D (pa - mf)[(p- k) 2 - m* 2 ] ' (16) 
Here we set trl to be 4 in the cutoff and 2 D I 2 in the 

s 

dimensional regularization schemes respectively. 

Substituting the solution of the gap equations (j4)) and 
(|5|) in the self-energy expression for 7r°, K° and K° 
and evaluating the pole structure of the denominator in 
Eq-®, we obtain the on-shell conditions for pion and 
kaon masses, 



A. 7r and K masses 

First, we consider pion and kaon. The masses of these 
mesons are obtained by observing the pole structure in 
their propagators. Employing the random-phase approx- 
imation (RPA) and the l/iV c expansion, the meson prop- 
agators are given by [3, [ 



A P (fc 2 ) 



2K n 



1 - 2K a l\ P {k 2 ) 



0(N C 



(9) 



= — £ - 2K 3 k 2 I uu (k 2 )\ k 2 =m 2 , 
m* 



= 1- 

-2G 



mz — m„ m,, — m„ 



2m* 2m* 
fitrS 11 - itrS s itrS s - itrS u 



V m l m *s 
2K 6 [k 2 - (m* s - m* u ) 2 ] J ns (fc 2 )| fc2=m 



(17) 



(18) 



These equations are used to determine the values of the 
constituent quark masses m* and m*. Because of the 
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SU(2) isospin symmetry the equations for the charged 
pion and kaon cases are the same as for the neutral ones. 
The QED effect is also important for the mass differences 
between the neutral and charged mesons. It should be 
noted that Eq. (fT8)) becomes Eq. (fl"7| in the limit m s — > 
m u because the SU (3) flavor symmetry is restored in this 
limit. 



B. n and K decay constants 

Next, we consider the pion and kaon decay constants. 
The decay constant fp is defined by the matrix element 
of the axial current between meson and vacuum states, 



= —M. 



4-D 



d D p 

{2n) D 



tr 



T 

luls-^S'gpqq^^TlS 3 



(19) 



where we introduce a renormalization scale Mq to de- 
fine the decay constant with a correct mass dimension, 
dim(/p) = 1 for dimensional scheme. The trace runs over 
flavor, spinor and color indices. The meson-to-quark- 
quark coupling gp qq is given by 



M, 



{ _ D dU P (k 2 



dk 2 



(20) 



Inserting Eqs. (JT5J| and (0 into Eq. ([20]). we obtain the 
pion and kaon decay constants, 



D-A 



1 



(21) 



Jus(0) 
r-1 

x / dx 



m*Jus(0) 
d D p 



trl • Nrlm*. 



ii 



where J us is defined by 



Ju S (k 2 ) =Iij{k z ) + 2N c {m* s - m*) 



* \2 



(22) 



dx 



d D p 



x(l — x) 



i(2n) D ( p 2-L us (k 2 )+ie) 3 - 



(23) 



Lij is defined in the appendix (see, Eq. (|A3jl ). As con- 
firmed in the case with m n and ttik, Eq.([22p corresponds 
to Eq. (f2~lj) in the m s — > m u limit due to the restoration 
of the SU(3) flavor symmetry. 



C. r\ and rj' mesons 

As is well known, the octet state, r\%^ and the singlet 
state, ?7q , are mixed in the real world. Thus the eigen- 
states, 77 and 77', with diagonal mass matrix are described 



as mixed states of 77s and 770- In the RPA the propagator 
of the 77 — 77' system is given by @, 0] 

A + (k 2 ) = 2K + [l-2K+Il(k 2 )}- 1 , (24) 

where K + and II are the 2x2 matrices 



K+ = 



n 



-^80 -^88 

rioo n 8 
ri8o n 88 



(25) 
(26) 



with 



Kqo = G — -K(i trS s + 2i trS u ), 

O 

K 88 = G--K(i trS s - 4t tiS u ), 
6 



Kq8 = K so 



V2 



K{itiS s -itrS u ), 



and 



noo(fc 2 )--[2nr(A: 2 ) + nr(fc 2 )] 
n 88 (fc 2 ) = |[nr(fc 2 ) + 2nr(fc 2 )] 



n 08 (fc 2 ) = n 80 (fc 2 ) 



2V2 



[nr(fc 2 )-nr(fc 2 )] 



To obtain the 77 and 77' meson masses we diagonalize the 
inverse propagator via an orthogonal transformation [3, 
Il4j |. Thus the rj and r( meson masses are given by the 
solution to the following equations, 



A(k 2 ) + C(k 2 ) - ^{A(k 2 ) - C(k 2 )} 2 + 4{B(k 2 )} 2 = 



for k = mz, 

(27) 

A{k 2 ) + C(k 2 ) + ^{A(k 2 ) - C(k 2 )} 2 + 4{B(k 2 )} 2 = 

for k 2 — m 2 , , 

(28) 

where 

A(k 2 ) = K 8S - 2II 00 (fc 2 ) det K+ , 
B(k 2 ) = -K QS -2Tl 08 (k 2 )detK + , 
C(k 2 ) = K Q0 - 2n 88 (fc 2 )detiv + . 

The mixing angle 9^ is found to be 

2B(fc 2 ) 



tan(20„) 



A(k 2 ) -C{k 2 )' 



(29) 



The off-diagonal matrix elements in Eq.(|24|) disappear 
and the mixing angle 9 V vanishes in the limit m s — > m u . 
In this limit Eg. ([27]) coincides with Eqs.(fT7| and (fig)) 
and the mass spectrum for the octet mesons degenerates. 
Since the Ua(1) anomaly breaks the degeneracy between 
the octet and the singlet mesons, a different on-shell con- 
dition is obtained from Eq. 
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D. Topological susceptibility 

To compare the QCD axial current with the NJL axial 
current we introduce the topological charge density Q(x) 



at 3, 4, 5 and 6MeV, and study the region, 2 < D < 4, in 
which an UV stable fixed point appears for G. Following 
the previous study Q, we fit the other 4 parameters {M , 
G, K, m s } by using the measured physical observables 
171, 



Q(x) 



32tt 



2 F^F afw = 2Klm[detq(l - l5 )q], (30) 



where g is the strong coupling constant of QCD and F" 
is the field strength for gluons. The topological suscep- 
tibility x is defined by the correlation function between 
the topological charge densities at different points, 



X = / d 4 x(0\TQ(x)Q(0)\0) 

connected- 



(31) 



It describes some global feature of QCD dynamics. In 
the leading order of 1/N C expansion it is given by [l8| 



4K 



X : 



(itrS u ){itTS s ) 



(2itrS s itrS u 
\ ml m* 



V6 



(2ztr5 s + ztr5 u )(n 00 (0),n 08 (0)) 



v/3 



(ttiS s - itr5 u )(n 08 (0),n 88 (0)) \ A+(0) 



±(2itvS s +itrS u )( 



itxS u ) 



n 08 (o) 
n 88 (o) 



(32) 



We compare this result with the one obtained by the 
lattice QCD. 



IV. PARAMETER SETTING IN THE 
DIMENSIONAL REGULARIZATION 

In this section we discuss the physical scale and the pa- 
rameter setting for the model with the dimensional regu- 
larization. In the study with the help of the NJL model, 
parameters are usually determined by fitting the physi- 
cal quantities {m^, /„■, mjr, X}, where there are several 
choices for X. Here we fit the parameters by choosing 
X = m v >. We reproduce these four observables in the 
dimensional regularization scheme without fixing two of 
model parameters, m u and D. As for the up quark mass, 
m u , we fix it by hand and test several values. The dimen- 
sion, D, is kept to be a free parameter. Then we calculate 
some meson characteristics as functions of D. Here we 
evaluate D-dependence of some observed physical quan- 
tities: the kaon decay constant fjc, r\ meson mass 
and topological susceptibility x- 



= 138MeV, m K = 495MeV, 
U = 92MeV, 77v = 958MeV. 



(33) 



These quantities are calculated by analyzing the Eqs. 
([17 )1 . ([15 ]) , PT|| and (J2SJ), respectively. We solve these 
equations under the constraints imposed by the gap equa- 
tions for m* and m* . Eliminating the term KitrS s from 
Eqs. (0} and (IT71) . we obtain m* as a function of m u and 
D, 



= 



1 771 

2~ itrS 



(34) 



We numerically evaluate this expression to plot m* as a 
function of D in Fig.[T] Inserting this m* , function of D, 
into Eq. (l2~Tj) . we also describe the renormalization scale 
Mq as a function of m u and D. The D-dependence of Mq 
is shown in Fig. [2] From Figs. [T] and[2l we confirm that 
m* and Mo in the 3-flavor NJL model have a behavior 
similar to that in the 2-flavor model [8j. 
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FIG. 1: The solution of the gap equation for m*i as a function 
of D for m u = 3, 4, 5 and 6MeV. 



Next we would like to fix the remaining 3 parameters 
{G, K, 777 s }. From Eqs. g|), © and JTSJ the coupling 
constants G and K are obtained as functions of m u ,777* 
and m* 



A. Parameter setting 

The model has 6 parameters m u , m s , G, K, Mq and D, 
as is already mentioned in Sec. ITT] Here we set m u (= nid) 



G{m* s ) 
K{m* s ) 



m* u ~m u 1 „ 
= Aim, 

9(m*) 



itrS s 



h{m*Y 



(35) 
(36) 
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FIG. 2: The renormalization scale Mo as a function of D for 
m u = 3, 4, 5 and 6MeV. 



where g(m*) and h(m*) are given by 



?K) = - i + 



2to? 



+ 



itrS 11 



itrS 8 - itrS y 



2itiS u 

+ {™?k ~ ( m *u " m s) 2 } Ius(m 2 K )] , 



itiS u itrS s 
— + 



h(m*) =(itrS s - itrS u ) 

ml ml 
+ {m 2 K - {m* u - m* s ) 2 } I us (m 2 K )} . 

Substituting Eqs. (05} and into Eq. jTH}, we de- 
scribe the current strange quark mass, m s , as a function 
of m„m* and m*. Then the rj' meson mass (f2"5| is also 
expressed by a function of ra„m* and m*. From the 
expression for m„/ and Eq. (|34[) we evaluate to* numeri- 
cally. 

The D-dependence of m* is shown in Fig. El We should 
note that qualitatively the solution for m* has a similar 
to to* behavior. Decreasing the dimension from four, we 
observe that the absolute value of to* goes down almost 
linearly. A different tendency takes place near the di- 
mension two, a discontinuity is observed at D ~ 2.5. A 
physical solution for to* does not appear, neither it is not 
seen for m* in Fig. [T] 

Once to* is obtained, G(m*), K{m* 8 ),m B {m*) are cal- 
culated by inserting the solution for to*. The numerical 
results for G, K and m s are shown in Figs. |H [5] and HI 
respectively. In Fig. 0] we observe a similar behavior of 
the four-fermion coupling, G, to one in the 2-flavor NJL 
model Q, while we do not have a consistent solution for 
the dimension where no solution is found for m*. There 
is no analog of Fig.[S]in the 2-flavor case. In Fig.[5Jwe see 
that K gets drastically large near the dimension where 
a consistent solution is lost. Finally in Fig. [5] it is found 
that the value of the current strange quark mass m s is 
roughly constant for the dimension larger than that cor- 
responding to the discontinuity point, but it falls down 
as D decreases below the discontinuity point. 

In Figs. |31 HIES and HI we observe a discontinuity around 
D ~ 2.5 where the behavior drastically changed. It 



comes from the fact that the values of the couplings G 
and K become divergent when m* s approaches m*. The 
self-energy 11" (fc 2 = m 2 ,) is also divergent at m 2 , = 
Ami 2 . Therefore m* has no physical solution around 
to* ~ m* ~ toy/2, as is numerically confirmed in Fig. 
[3] This is the reason why the behavior of to* is different 
from that of to* near the dimension two. The constituent 
strange quark mass, to*, can not develop a value smaller 
than to ?? '/2. At the discontinuity the rj' meson propaga- 
tor contains a imaginary part which corresponds to the 
decay width. 
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FIG. 3: The solution of the gap equation m* u as a function of 
D for m u = 3, 4, 5 and 6MeV. 
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FIG. 4: The 4-fermion coupling G as functions of D for m u 



B. Physical quantities 

We have fixed the model parameters in the previous 
subsection except for the up quark mass m u and the di- 
mension, D. We are now ready for evaluating various me- 
son properties as a function of D for a fixed m u . Employ- 
ing the obtained parameters and the chiral condensates, 
to* and to* , one can numerically calculate the kaon decay 
constant, fx, the r\ meson mass, to, ( , and the topologi- 
cal susceptibility, \. The behavior of fx, Tn„ and \ are 
displayed in Figs. [3 [5] and [SI respectively. These values 
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FIG. 5: The 6-fermion coupling K as functions of D for m u — 
3, 4, 5 and 6MeV. 
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FIG. 6: The strange quark mass m s as a function of D for 
m u = 3,4,5 and 6MeV. 
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FIG. 8: The r\ mass m v as a function of D for m u — 3, 4, 5 
and 6MeV. 
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FIG. 9: The topological susceptibility \ as a function of D 



are roughly constant near four dimensions. They rapidly 
fall down as D decreases near two dimensions. We again 
observe a discontinuity around m* ~ m* ~ The 
topological susceptibility blows up around the disconti- 
nuity. 
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FIG. 7: The pion decay constant fx as a function of D for 
m u = 3, 4, 5 and 6MeV. 

As is seen in Fig. [7J the decay constant is smaller 
than the observed one, HOMeV, in the region, 2 < D < 
4. The model has to be improved to describe the Kaon 



decay. We can fit the r\ meson mass to m n = 548MeV 
by tuning the dimension of the fermion loop integrals, 
D. The topological susceptibility is calculated to be 
X 1/4 = 170±7, 174±7MeV in lattice simulations [H| and 
y 1 / 4 = 179MeV in the Witten-Veneziano mass formula 
[20l [2l1 ] . We can also fit the topological susceptibility in 
the interval, 2 < D < 4. Consequently the dimension for 
the fermion loop integrals, D, can be fixed to reproduce 
the 77 meson mass or the topological susceptibility. We 
will discuss this matter in the next section where we need 
to fit an additional parameter D. 



V. QUARK MASS DEPENDENCE 

In the cutoff regularization all the model param- 
eters can be fixed by fitting the physical quantities 
ttik, "vl for a given m u . In the dimensional 
regularization an an additional observable, X' , needs to 
be fitted due to one more parameter involved in this case. 
Therefore, there are five physical quantities, i.e. the set 
of them, {m^, f„, mx, ""V, X'}, which are needed to 
fit all the model parameters. In this section we con- 
sider two different observables for X'; one is the rj meson 



mass, X' 



another is the topological susceptibil- 
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ity, X' — x- By using the experimental/empirical values 
of m n and x we can define two different parameter sets 
and evaluate the meson properties. Then some physical 
quantities are evaluated for m u — 3, 4, 5, 5.5 and 6MeV. 
We calculate the physical quantities in the both regular- 
izations to compare them with each other for each given 



A. Physical quantities in the dimensional 
regularization 

To show the validity of the model as a low en- 
ergy effective theory of QCD here we evaluate pa- 
rameters of the model to reproduce the input sets 
{m W) U, m K , m v >,m v } or {m^, /„., m K , m v >,x} in the 
dimensional regularization and discuss the other (out- 
put) physical quantities. There exist two solutions to 
reproduce m n = 548MeV for m u < 5MeV. In Table U 
we show m s , fa, X, (uu) r under the condition that 
rn v = 548MeV. The subscript r in (uu) r stands for 
that the quantity is renormalized, (uu) r = (uu)Mq~ d . 
The last lines in Table Q] and |TT] show the experimen- 
tal/empirical values. The light quark masses are eval- 
uated at lGeV [13] ■ The empirical value of (uu) ~ (ss) 
is evaluated by using the Gell-Mann-Oakes-Renner rela- 
tion (22|. We observe that the topological susceptibility 
X is always larger than the one in the lattice simulation 
and the Witten-Veneziano mass formula. 

In the Tables |TT] and IIIII we fix the parameters of the 
model so that to reproduce x 1 ^ 4 — 170 or 179MeV 
and then we calculate m s , Jk, m v , (uu) r for m u — 
3, 4, 5, 5.5, 6MeV. As is shown in Fig. |H1 we find two 
solutions in both sides of the discontinuity to satisfy 
X 1 / 4 = 179MeV for m u < 4MeV. In any case, we obtain 
the value for the r\ meson mass which is smaller than the 
observed one, m n = 548McV. 



TABLE I: Physical quantities in the dimensional regulariza- 
tion (in units of MeV, except for D), is fixed at 548MeV. 
The last line shows the experimental/empirical values. 



m u 


m s 


Ik 


x 1/4 




D 


3.0 


84.9 


90.8 


244 


301 


2.289 


3.0 


79.0 


91.4 


225 


301 


2.372 


4.0 


118 


91.4 


254 


274 


2.378 


4.0 


106 


92.1 


225 


273 


2.522 


5.0 


156 


92.0 


275 


254 


2.465 


5.0 


134 


92.7 


224 


254 


2.687 


5.5 


148 


92.9 


224 


246 


2.775 


6.0 


162 


93.2 


224 


239 


2.867 


3.4 - 6.8 


94.5 - 176 


110 


170 - 179 


228 - 287 





TABLE II: Physical quantities in the dimensional regulariza- 
tion (in units of MeV, except for D), x 1/4 is fixed at 170MeV. 
The last line shows the experimental/empirical values. 



m u 


m s 


Sk 


m, 


/- \l/3 
— [UUjr 


D 


3.0 


77.1 


88.4 


481 


302 


2.280 


4.0 


105 


88.8 


478 


274 


2.360 


5.0 


134 


89.2 


475 


255 


2.434 


5.5 


150 


89.3 


473 


247 


2.467 


6.0 


166 


89.5 


471 


240 


2.500 


3.4-6.8 


94.5 - 176 


110 


548 


228 - 287 





TABLE III: Physical quantities for the dimensional regular- 
ization (in units of MeV, except for D), y}^ is fixed at 
179MeV. 



m u 


m s 


}k 


m v 


— {UU) r 


D 


3.0 


78.0 


88.7 


498 


301 


2.281 


3.0 


74.8 


92.7 


507 


300 


3.222 


4.0 


106 


89.1 


495 


274 


2.363 


4.0 


100 


92.9 


507 


272 


3.895 


5.0 


136 


89.5 


491 


255 


2.438 


5.5 


152 


89.7 


489 


247 


2.472 


6.0 


168 


89.9 


487 


240 


2.505 


B. 


Physical quantities 


in the cutoff regularization 



In this section we fix the model parameters to repro- 
duce the measured values of Eq. ([33)) in the cutoff regu- 
larization for a given m u . The quantities, {m s , G, K, A}, 
are found following the same procedure as in the previ- 
ous section. In the Table IIVI we show the parameters, 
G,K and A for m u = 3, 4, 5, 5.5 and 5.87MeV. It is ob- 
served that the cutoff scale A decreases as m u increases. 
No solution is found to simultaneously satisfy Eqs. (fTT|) 
and (|2l|) for m u > 5.87MeV. In Table [V] the quantities, 
m ril Xi m s, (uu), are shown. In the cutoff regulariza- 
tion the kaon decay constant, fx, is consistent with the 
observed value for a smaller m u , while the rj meson mass 
and the topological susceptibility are smaller than the 
experimental/empirical values. We obtain only smaller 
rj meson mass than the experimental one, similarly to 
the situation in the dimensional regularization with fixed 
topological susceptibility x- 

Thus none of the regularizations can describe all of the 
nonet meson properties. We can fit the rj meson mass or 
the topological susceptibility in the dimensional regular- 
ization but we obtain only a smaller kaon decay constant 
in either case. In the Pauli-Villars regularization it has 
been found that m u — rrid = 2.7MeV, m s — 92MeV, 
fjc = 131MeV and m v = 526MeV where the input pa- 
rameters rri-x , rrtK, and rrirf have been used [231 ]. 
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TABLE IV: Parameter list for cutoff scale A, 4-fermion cou- 
pling G and 6-fermion coupling K. 



m u (MeV) 


A(MeV) 




GA 2 


KA 5 


3.0 


960 




1.55 


8.34 


4.0 


797 




1.60 


8.38 


5.0 


682 




1.71 


8.77 


5.5 


630 




1.81 


9.17 


5.87 


580 




2.09 


10.1 


TABLE V: Physical quantities 
units of MeV). 


in the cutoff regularization (in 


m u m s 


Ik 


m n 


x 1/4 




3.0 89.5 


113 


451 


160 


301 


4.0 110 


107 


457 


158 


273 


5.0 128 


101 


473 


160 


253 


5.5 136 


97.3 


482 


163 


245 


5.87 139 


93.3 


501 


172 


240 



other hand, the r\ meson acquires a smaller mass than the 
observed one for the fixed x- 

In the cutoff regularization the cutoff scale significantly 
depends on m u ; the cutoff A increases with decreasing 
m u . It is interesting to note that the coupling constants 
become smaller when one takes the larger cutoff, which 
is consistent with the renormalization group argument 
where the coupling strength becomes smaller with in- 
creasing the energy-scale. This cutoff effect is also con- 
firmed numerically in the NJL model through changing 
the cutoff in the temporal direction [24| . 

In the Tables. [Til Mil and [V] some difficulty is seen 
in tuning Jk, and x simultaneously. Thus, this pa- 
per shows that the model based on the Lagrangian ([T]) 
is not satisfactory in either of regularizations. It teaches 
us that new terms should be added to the Lagrangian 
to improve the model at least in one of the regulariza- 
tions. Especially, we are interested in including vector 
type and multi-fermion interactions and hope to report 
on the problem in future. 



VI. CONCLUSION 

We studied nonet meson properties in the three-flavor 
NJL model with the dimensional and sharp cutoff reg- 
ularizations in the leading order of the 1/N C expansion. 
We employed m u , m^, itik, f^, m n i as input parame- 
ters and fix the model parameters, m s , G, K, Mq, and 
m s , G, A in the dimensional and cutoff regulariza- 
tions, respectively. In the case of the dimensional reg- 
ularization the dimension, Z?, is still a free parameter. 
Thus we evaluate the kaon decay constant, the r\ meson 
mass and the topological susceptibility as a function of 
D. 

The constituent up-quark mass, to*, and the renor- 
malization scale Mq behave in a similar to the two-flavor 
case way. No consistent solution for m* is found and 
the effective coupling, G and K, are divergent around 
to* ~ to* ~ toy/2. Below this region rj' meson propa- 
gator develops an imaginary part. The kaon decay con- 
stant, /if, is smaller than the observed value, HOMeV, 
for 2 < D < 4. 

We found that m s , /k, m n ^ x are almost constant 
near four dimensions. Fitting the parameters at lower 
dimension and then taking the four dimensional limit, ob- 
tained values for m s , fa, rn rn x are not divergent. The 
results do not depend on the regularization parameter. 
It bring us an interesting idea of making a regularization 
independent prediction in the NJL model. 

We have evaluated the physical quantities, m s , fx, X, 
m v and (uu) , in both the dimensional and the cutoff reg- 
ularizations for a fixed m u . In the dimensional regular- 
ization the topological susceptibility, x, develops a larger 
value than the one obtained in the lattice simulation and 
Witten-Veneziano mass formula for the fixed m n . On the 
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Appendix A: Integrals 

In this section, we demonstrate several integrals both 
in the cutoff and dimensional regularization schemes 
which are required in the section [TT] and IIIII 



1. The cutoff regularization 

The chiral condensates i trS 1 in Eq. © are important 
quantities because they concern gap equations. They 
take the following form after the integration is performed, 

= ^<{ a V /a2 + < 2 

, 2 A + y/A 2 + TO* 2 | 

V m i J 

The next integral to carry out is Iij(k 2 ) in Eq. tfll))) which 
determines self-energies of mesons, lip, Eq. (IT2l . It takes 
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the form 



In = 



where 



Nc 
8tt 2 J 



dx { —- 



A 



^A 2 + L lJ (fc 2 



In 



A+ y/A 2 +i, tJ (fc 2 ) 



(A2) 



Lij{k 2 ) = to* 2 - (mf - mf)x - k 2 x{l - x). (A3) 

Finally, Eqs. (j2"2")l and ([25)1 involved in the derivation of 
the kaon decay constant are calculated as 



Regarding i^, we have different expressions depending 
on the value of fc 2 ; 



Nr. 



(4tt) d / 2 



D/2-2,,2 



(fc^), (A9) 



for (m* - to*) 2 < k 2 < (to* + to*) 2 and fc 2 = 0, 



J«.(0) 

i- 1 



m* u I us {0) - ^(m* - m*) 



A 



o 



v/A 2 + L us (0) 



In 



A +^ + ^(0) 



V^(0) 



Jus{k ) — Ius{k ) 



AT C 

167T- 



rix 



.a 3 (to:-t<) 2 

a;(l — x) 



o L us (k 2 )[K 2 + L us {k 2 )Y/ 2 



(A4) 



(A5) 



Zy(fc 2 ) 



A T , 



x < a 



D/2-1 



D/1-2 



D D 

F 2 , 1, 

2 ' 2 '2 



13-2 ,J 



|1 -" + )" L> F ( 2 -f'f 1 



L> 1 - a. 



2 ' a+ — a_ 



+ e ^-D/2) {k 2 ) l-D/2 l/ D-3 B _ _ \ 



(A10) 



2. The dimensional regularization 

The corresponding integrals in the dimensional regu- 
larization are performed as 



for < k 2 < (to* — to*) 2 , (to* + to*) 2 < fc 2 , respectively. 
F denotes the hypergeometric function. In Eq. (|A10I) we 
introduce the quantities a± and i/y which are defined by 



itrS 1 

S,C 



Nc r fi- — to^to* 2 ) 15 / 2 1 

(2^)^/ 2 V 2 ' a 1 J 



(A6) 



2 D/2 



M ^- 4 J„ s (0) 



m:/ Ufl (0)- 



(4tt)^/ 2 I 2 



x(m:-m*) / dxxL us (Q) D/2 ~ 2 
Jo 



(A7) 



N c ( D\ , 
rTT^r 3- — (m s -m^ 



(4tt) d / 



x / <ia;a;(l — x)L ua (k ) 
Jo 



2\D/2-3 



(A8) 



a± = 



fc 2 + (to* 2 - to* 2 ) ± 



2fc 2 



Vij = - 2fc 2 (m* 2 + to* 2 ) + {mf - to* 2 ) 2 
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